The Molecular Universe 

Proceedings I A U Symposium No. 280, 2011 

Jose Cernicharo and Rafael Bachiller, eds. 



© 2011 International Astronomical Union 
DOI: 00.0000/XOOOOOOOOOOOOOOOX 



Radiative transfer and molecular data for 

astrochemistry 

Floris van der Tak 

SRON Netherlands Institute for Space Research, Landleven 12, 9747 AD Groningen, The 

Netherlands 
email: vdtak@sron.nl 

l" Abstract. The estimation of molecular abundances in interstellar clouds from spectroscopic ob- 

,_^ ' servations requires radiative transfer calculations, which depend on basic molecular input data. 

C^ , This paper reviews recent developments in the fields of molecular data and radiative transfer. 

^N) ' The first part is an overview of radiative transfer techniques, along with a "road map" showing 

I which technique should be used in which situation. The second part is a review of measurements 

^ , and calculations of molecular spectroscopic and coUisional data, with a summary of recent coUi- 

' • sional calculations and suggested modeling strategies if collision data are unavailable. The paper 

^\ concludes with an overview of future developments and needs in the areas of radiative transfer 

and molecular data. 
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1. Introduction 

At the heart of observational astrochemistry is the estimation of molecular abundances 
from telescope data. Extracting the information out of the data requires radiative transfer 
calculations which need molecular input data, and these two topics are discussed in this 
paper. In many cases, measuring the chemical composition of interstellar gas is coupled 
to the estimation of its physical conditions, in particular the kinetic temperature and 
the volume density. In addition, molecular abundances usually need interpretation with 
a chemical model, to estimate parameters such as the age of the cloud or its ionization 
rate by cosmic rays. These two important related topics are not covered here. 

Before starting the review, it is useful to consider the limitations of radiative transfer 
calculations. Some limitations are of observational nature, such as resolution mismatch 
between observations at different wavelengths and the presence of unresolved substruc- 
ture in the data. In addition, the lines in a spectrum may not all sample the same gas, 
for instance when some lines appear in emission and others in absorption. Furthermore, 
the observations may only cover part of the source, and only a few molecular energy 
levels. Finally the radiative transfer calculations may be limited by the availability of 
molecular data: either spectroscopic information or collisional excitation probabilities. 
When interpreting observations, it is important to keep these limitations in mind and to 
choose the model such to minimize their impact. 

The field of molecular data and radiative transfer is large, and due to space limita- 
tions, this paper ca nnot treat every aspec t in depth. This review focuses on the analysis 
of spect ral lines; see Pascucci et al.l (J2004I ) for a discussion of continuum radiative trans- 
fer and IShirlev et al.l (120111) for a recent update about dust opa cities. In the area of 



line ra diative transfer, see iPeraiahl (|200ll) for an introduction and iPavlvuchenkov et al 



(120071) for re c ent d evelo pments. For detaile d accounts on laboratory spectrosc opy, see 
Pickett et al.l ( 1998 ) and lMiiller et al.l ( 2005 ) and references therein. The book by lFlower 
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( 20031 ) reviews the field of interstellar molecular colli sions. Finally , the estimation of phys 



ical parameters of interstellar clouds is reviewed by Evand (|l999i ) . Note that throughout 



this review, the term 'molecules' includes mono-atomic and ionic species. 

2. Analyzing molecular line spectra 

The goal of radiative transfer calculations is to estimate physical parameters from 
astronomical observations. In the case of astrochemistry, the prime application of radia- 
tive transfer is the estimation of molecular abundances from line spectra at infrared and 
(sub)millimeter wavelengths. This review focuses on the analysis of radio and far-infrared 
data, but most of the mentioned techniques are also applicable into the mid-infrared. The 
near-infrared, visible and ultraviolet spectra of circumstellar envelopes and disks are also 
of interest in astrochemistry and make additional demands on molecular data and radia- 
tive transfer; however, these issues are beyond the scope of this paper. 

There is no single best method for radiative transfer calculations; rather, the amount 
and nature of available data and other constraints determine the appropriate level of 
analysis. For example, the use of non-LTE models is often limited by the availability 
of molecular collision data (^ 13. 2p . Also, applying a two-dimensional model with 1000 
grid points is clearly overkill if only a handful of spectral lines have been observed in a 
single direction. On the other hand, the increasing data volumes from modern telescope 
facilities demand ever faster and more sophisticated models for the data analysis, and 
developing such tools is an active area of work. This section reviews radiative transfer 
models at three different levels, ranging from fairly basic LTE models to fairly complex 
non-local models in 2D or even 3D. 

2.1. Single excitation temperature ('LTE') models 

The simplest way to estimate a molecular column density from an observed line intensity 
is to avoid radiative transfer altogether, and the first step in doing so is the estimation 
of the molecular column density in the initial state of the observed line (upper state 
for emission lines, lower state for absorption lines). The conversion from line intensity to 
optical depth is straightforward for absorption lines, but involves assuming a beam filling 
factor if the lines appear in emission. Observations at higher angular resolution help to 
constrain this factor if they are available; otherwise a point source and a very extended 
source provide useful limiting cases. 

The second step is to account for the other, unobserved states of the molecule, and 
the common approach is to assume that the energy levels are populated according to 
the Boltzmann distribution. The characteristic 'excitation temperature' (Tox) is a priori 
unknown, but is often taken to equal Tkin, which is the case of local thermodynamic equi- 
librium (LTE) which applies at high densities when collisions determine the excitation. 
At low densities, radiative decay competes with collisions, so that generally Tox < ^kin, 
but even in this case, this method is often (misleadingly) called the LTE method. 

If several lines from a range of energy levels of the same molecule have been observed, 
their relative intensities can be used to constrain Tex- The most direct way to do so is to 
fit a straight line to the line intensities as a function of upper level energy. When rota- 
tional lines are probed, the popular name 'rotation diagram' is correct, but in general, 
the alternative name 'Boltzmann plot' is more appropriate. For low line optical depths, 
this method yields reliable est imates of Tpx c^nd mole cular column density, both for sub- 



millimeter emission lin es (e.g., Kristensen et al.ll2010l ) and for visual/infrared absorption 



lines fe.g.. lKazmierczak et al.ll2010 ). However, at high column densities, the derived T( 



is inaccurate because of optical depth effects, and correction factors must be applied 
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( Goldsmith fc Langeiill999l ). The advantage of this improved method (sometimes caUed 



the 'population diagram') is that it also provides an estimate of the size of the emit- 
ting region. An example of combining the two methods is the analysis of CH3OH line 
images of the Cep A region, where classic rotation diagrams work well for the diffuse 
extended gas, while opt ical depth correction are needed close to the central protostar 



( Torstensson et al.ll2011 ) 



Many programs to perform 'LTE' calculations exist, but one popular program is CAS- 
SIS which is publicly availablqjj. Besides multi-line models for single species, this program 
is capable of synthesizing wide-band spectra based on models for many molecules, each 
with their own set of excitation conditions. 

2.2. Non-LTE models 

The LTE method is fast and applies to any molecule with a known spectrum, but its 
assumption of a Boltzmann distribution for the level populations often does not hold, 
and measuring Tox is often difficult. More sophisticated non-LTE methods retain the 
assumption of a local excitation, but solve explicitly for the balance between excitation 
and de-excitation of molecular energy levels. Both collisional and radiative processes 
contribute to this balance, and these 'non-LTE' models require molecular collision data 
as input, in addition to spectroscopic data. This requirement limits the use of non-LTE 
models as collision data do not exist for all astrophysically relevant species. However, 
the advantage of non-LTE models is that not only the column density but also the 
kinetic temperature and the volume density of the gas can be constrained from multi- 
line observations. 

There are two main ways in which non-LTE models treat the effect of nonzero line 
optical depth. One way is the Sobolev approximation, where the cloud is assumed to 
have a significant velocity gradient, of the order of 1 kms~^ pc~^ or more. The resulting 
Dopplcr shift of the line profile decouples the excitation at one end of the cloud from radi- 
ation from the other, and has the effect of localizing the radiative transfer problem. Such 
models match observed spectra, but the reliability of the derived p arameters is unc lear 
since most clouds do not have systematic velocity gradients (e.g.. IOssenkopall997f ). A 



somewhat better treatment of optical depth effects is the 'escape probability' formalism, 
in which the chance for photons to leave the cloud is based on the current local estimate 
for the optical depth and an assumed geometry. Both treatments require iteration be- 
cause the molecular excitation and the radiation intensity are mutually dependent, but 
on modern computers, convergence is fast. 

Many programs to perform non -LTE calculations exist, but one popular program is 
RADEX ( Van der Tak et al.l 12007 ) which is publicly available^. This program exists as 



an on-line calculator for quick estimates of line intensities from assumed column densities 
or vice versa, as a module inside the CASSIS package, and as a stand-alone program for 
extensive calculations, for instance over a grid of kinetic temperatures an d volume den- 
sities. An example of this latter use is the work bv lCinsburg et al.l ( 2011 ) who constrain 



the densities of molecular gas close to ultracompact H II regions with non-LTE models 
of H2CO line ratios. 

Multi-zone non-LTE models divide a cloud into grid cells and calculate the escape prob- 
abilities for photons from one cell into each of the others. These programs are much faster 
than true non-local radiative tr ansfer programs but are often limit ed to plane-parallel 
geometry ( Elitzur fc Asensio Ramos 2006) or particular molecules ( Poelman fc SpaansI 



t |http: //cassis. cesr.fr/ 1 

J |http : //www ■ sron . rug . nl/ ~vdt ak/r adex/index . shtml | 
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20051 ). Versions of the m ulti -zone escape probability m ethod in spherical geometry exist 
as well (JMalonevI llQSTt see Ivan Dishoeck et al.l Il99ll). and the re sults appear to be in 



good agreement with those of non-local methods ( Yun et al.ll2009l ) 



2.3. Non-local models 

The most sophisticated radiative transfer programs use non-local methods: they drop the 
assumption of a local excitation and solve for the molecular excitation as a function of 
position in the cloud. To estimate the local radiation field at all line frequencies due to 
radiation from everywhere else in the cloud, two main methods exist. Accelerated Lamda 
Iteration (ALI) program discretize the cloud onto a grid of points and split off the local 
contribution from the full radiation field at each point. Monte Carlo (MC) methods 
divide the cloud into grid cells and send photon packages in random directions from each 
cell. The performan c e and convergence of bot h types of programs has been tested by 
Van Zadelhoff et al.l (|2002[ ) for HCO+ and by IVan der Tak et al.l (|2005l) for H2O. The 
high accuracy of these models comes at a price: for complicated model geometries, large 
gradients in excitation conditions, and molecules with many energy levels, the calculations 
can take several hours. 

Several programs to perfor m non-local radiative tr ansfer calculations exist, but one 
popular program is RATRAN (JHogerheiide fc van der Tak 2000 ). The ID version of this 
Monte Carlo program is publicly availabl qil and the 2D version is availa ble upon request. 
Recently, the program LIME has appeared p Brinch fc Hogerheiidd uOloV which is capable 
of doing 3D non-local calculations. This program is not publicly available, but the authors 
provide a copy upon requesly. 

Programs such as RATRAN and LIME are useful to constrain the physical (tem- 
perature, density and velocity) structure of astrophysical sources, and also to estimate 
molecular a bundances as a fun ction of position in a cloud. An example of the former is 
the work bv lPanic et al.l ( 20091 ) . who model interferometric images of IM Lup in the CO 
J=2-l line to constrain th e 2D density stru c ture o f its protoplanetary disk. An example 
of the latter is the work bv lChavarria et al.l ( 2010l ). who model multi-line Herschel-HIFI 
observations of II2O toward the high-mass protostar W3 IRS5 to constrain the distribu- 
tion of H2O within the protostellar envelope. 



3. Molecular input data 

The interpretation of interstellar line spectra with radiative transfer calculations gen- 
erally requires two kinds of molecular input data: spectroscopic data (energy levels, sta- 
tistical weights, line frequencies, and Einstein A coefficients) and collision data (for all 
but the simplest LTE models). This section describes how such data are obtained, what 
the results are, and which limitations exist. 

3.1. Spectroscopic data 

Characterizing the spectrum of a molecule begins with measuring its line frequencies in 
the laboratory. Many species are available commercially; many others can be produced 
in the reaction between two or more species. The molecules are let into a gas cell, which 
is typically a few meters long and made of glass. The gas cell is pumped to high vacuum, 
a radiation source is put at one end and a spectrometer at the other. 

The bandwidth over which the line frequencies are recorded depends on the measure- 
ment technique, and is limited by atmospheric transmission. Ideally, spectra are measured 



f 'http : //www . sron . rug . nl/~vdtak/ratraii/ 
:j: http: //www. strw. leidenuniv.nl/~brinch/website/limecode .html 
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in the microwave, (sub)niilliineter, and near- and mid-infrared ranges. The resolution is 
usually limited by Doppler broadening of the lines, which amounts to i? = Aiy/i/o ^ 10^^ 
or 0.1 MHz at 100 GHz for air at 300 K. Special techniques to achieve sub-Doppler reso- 
lution exist but are not always available. 

The second step is to apply a Hamiltonian model and to assign quantum states to 
the measured lines. Even for linear molecules, this step is not trivial, since absorption 
spectra at room temperature usually contain several vibrational states, and sometimes 
multiple isotopologues may contribute. For large molecules with many vibra t ional modes, 
the number of parameters in the Hamiltonian becomes large; see lXu et al.l (J2004( ) for an 
example with internal rotation with wSO parameters. 

The standard program to perform Hamiltonian model fits is spf it by Herb Pickett; 
for certain specialized situations, other programs exist. Besides line frequencies and as- 
signments, such programs need the molecular dipole moment, which is oft en calculated 
ab initio, or sometimes measured directly fe.g.- lBotschwina fc OswaldlbOOSI ) . The output 
from these programs is a list of the complete molecular line spectrum, indicating where 
this prediction deviates from the measurements. These deviations may then be used to 
change some assignments in order to optimize the overall quality of the fit; 

The results of the spectroscopic work are collected in databases: JPl|j] (jPickett et al 



1998[ ) and CDM^ (iMiill er et ai]l2005l) are th e leading efforts at submillimeter wave- 
lengths and HITRAlNljl ( Rothman et al.l 120091 ) at infrared wavelengths. As of Summer 
2011, the submillimeter databases contain well over 500 entries, including many species, 
isotopologues, and vibrational states. The HITRAN database covers fewer species than 
JPL and CDMS, but many more vibrational states, as this database is aimed at atmo- 
spheric research. 

Although the completeness and accuracy of the current spectroscopic work is sufficient 
for many astrophysical applications, areas exist where astrochemical work is limited by 
spectroscopy. One example are the many unidentified fe atures in spectral line surveys, 
especially at frequencies >1 THz fe.g.. lBergin et al.ll2010l) . In this regime, line frequencies 
cannot be measured directly in the lab, so that their identification relies on Hamiltonian 
predictions which are sometimes quite uncertain. 

For large molecules with significant internal structure, a Hamiltonian fit to the labo- 
ratory spectra and assigning each of the man y lines from the variou s bending / torsional 
/ vibrational states is impractical. Instead, iFortman et al.l ( 20101) have measured the 
spectra of several such molecules at various temperatures so that features in spectral line 
surveys can be assigned to particular species, if not to specific transitions. Such efforts are 
essential for t he reliable identifica tion of new molecular species in crowded astronomical 
spectra (e.g.. lBelloche et al.ll2008l ). 

In some cases, astronomical observations are useful to characterize molecular spectra. 
The intrinsic line width of cold, quiescent interstellar cloud cores is smaller than achiev- 
able in the laboratory, so that higher spectral resolution can be obtained. In addition, 
certain species are difficult to produce on Earth or very sho rt-lived under terrest rial con- 
ditions. Examples include the sp ectra of N2H"'" and NzD"*" (jPagani et al.ll2009al ). DCO"*" 
(CaseUi fc Pore 2005) and DNC (|Van der Tak et al.ll2009l) . 



3.2. Collision data 

Besides spectroscopic data, modeling interstellar line spectra requires knowledge of the 
probabilities for (de-)excitation of the molecule upon collision with other molecules. In the 



f |http : //spec .j pi. iiasa. gov/ftp/pub/catalog/c atform .html I 
I http://www. astro . uni-koeln . de/cdms7 
^ http : //www . cf a . harvard . edu/hitran/ 
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case of interstellar clouds, the most common collision partner is H2, but collisions with H, 
He, H+ and electrons are also sometimes relevant. Collisional excitation probabilities are 
usually cast in the form of rate coefficients (with dimension cm^ s~^) which are collisional 
cross-sections (cm^) integrated over velocity (cms~^). These collision data are the results 
of quantum- mechanical calculations, with the occasional experiment for verification. Only 
downward collisions need to be calculated: the upward rates follow from detailed balance. 

The calculation of molecular collision data is a two-step process: first the potential 
energy surface of the collisional system is calculated, and second the dynamics of the 
interaction are evaluated, taking the Maxwell average of the velocity distribution. These 
calculations are time-consuming, which is why early work (before Ri2000) concentrated 
on closed-shell molecules, used He as collision partner, and adopted simplified methods 
for the calculation of the potential energy surface and the dynamics. Thanks to recent 
advances in computer power and generous funding from the European Union, modern 
calculations usually consider H2 as collision partner, use accurate potential energy sur- 
faces and close-coupling level dynamics, while treatment of open-shell molecules such as 
CN is now also possible. 

Experimental verification of collisional calculations is sometimes possible, albeit for 
a limited number of states. Besides pressure boadening experiments (see § 13.41) which 
measure the quality of the potential energy surface as a whole, detailed state-to-state 
experiments ar e valuable tests of qu antum- mechanical computations. One recent example 
is the study of lAsvanv et al.l ( 20091 ) on the H2D+ - H2 system, where both inelastic and 
reactive collisions play a role. 

Due to the computational limitations mentioned above, the field of molecular collisions 
is not nearby as adva nced as that of spe ctroscopy. As of Su mmer 2011, the two le ading 
databases, LAMDAQ ^Schoier et al.ll2005i ) and BASECOlU (|Dubernet et al.ll2006b[ ) con- 
tain collision data for 32 molecules. This number is only a fraction of the «160 molecules 
species k nown to exisg in int erstellar and circumstellar media (H2O2 being the latest 
addition; iBergman et al.ll201ll) . Since the detection rate of interstellar molecules is simi- 
lar to the calculation rate of collision data (a few molecules per year), this gap is unlikely 
to close - in fact it may even widen with time. Moreover, coverage of isotopologues and 
vibrational states is usually limited to the most common species in the ground state. 



3.3. Recent collision calculations 



Since the previous review on this topic (jVan der Tak fc Hogerheiidd 120071 ). many new 
calculations of molecular collision cross sections have appeared in the literature. Below 
is a list of calculations from the past 3 yea rs, which may not be complete. 

The calculations for H2D+ and D2H+ bv lHugo et al.l ( 20091 ) are a major step forward, 
even if they are based on a statistical approach rather than quantum scattering calcu- 
latio ns. These ions are crucial for th e study of cold and dense pre-stellar cloud cores 
(^Caselli et al.ll2008l : IParise et al.ll2011 ). Some of the state-to-state rates have been veri- 
fied by e xperiment, which is important for the study of the ortho-para ratio of H2, for 
example (IPa gani et al.ll2009b ). 

New collision rates for H2CO with H2 have been published bv lTroscompt et al.l ( 2009I ). 
which are valuable because H2CO line ratios are often used to measure the kinetic temper - 
ature (and sometimes the volume density) of interstellar clouds fe.g.. lMuhle et al.ll2007r ). 
The new rates are however incomplete as only the ortho species is covered; hopefully the 
rates for para-H2C0 will follow soon. 



f |http : //www ■ str w . leidenuniv . nl/~moldata/ 1 

I http://basecol . obspm . f r/ 

^ http: //www. astro .iini -koeln.de/cdms/molecules 
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Th e calculation of collision rates of HCN and HNC with He ( Sarrasin et al.lbOlOtlDumouchel et al 
20101 ) are useful because the HCN/HNC ratio is often used as a chemical diagnostic 
of temperature ( Loenen et al.ll2008t [Roberts et al.l 12011 ). Previously, HCN rates were 



used for HNC, which turns out not to be accurate at all. Recent observations of the 
H^'^CN/HN-'^'^C ratio in a sample of dense cores are in good agreement with the newly 
calculated collision rates ( Padovani et alj 2011b ). 

The collision rates for O2 with He bv lLigud ( 20101 ) have appeared ju st in time for their 

use to interpret the definitive detection of interstellar O2 with Herschel ( Goldsmith fc the HOP teami 
20111 ). This result is an im portant confirmation of the earlier tentative detection with 
Odin (JLarsson et al.l 120071 ). Although the excitation of O2 is usually not far from LTE, 
the new rates are useful to quantify this effect. 

New results have also appeared for the CN radical (ILigue fc Klos 2011 ), which is ofte n 
used as a probe of dissociating ultraviolet radiation (e.g.. IPerez-Bcaupuits et al.ll2009l ). 
especially together with HCN. The new collision rates are especially useful because they 
include the hyperflne structure of CN, which allows quantitative estimates of the optical 
depths of its emission lines. 

Models fo r high- J lines of CO have become much more accurate with the collision cal- 
culations bv lYang et al.l ( 20101) . These data are important for example to model Herschel 
observations of the CO emission from young stellar o bjects which heat their su rroundings 
to several lOOK by UV irradiation and shocks (e.g.. lVan Kempen et al.ll2010l ). 

One key interstellar organic molecule is CH3OH , and its lines are a diagnostic of 



both the kinetic temperature and the vo lume density of interste llar clouds (jLeurini et al. 



2007t ). The new collision calculations bv lRabli fc Floweij (|2010[) bring models of CH3OH 
excitation to higher accuracy. The coverage of the first torsionally excited state also 
enables to constrain the intensity of the far-infrared radiation field within astrophysical 
sources. Torsionally exc ited emission has already been observed in some objects (e.g., 
Torstensson et al.ll201ir) and will become routine with ALMA. 

The excitation of interstellar H2O is a long-standing problem in astrophysics. While 
the main iso tope is treated in the foll owing section, w e note here the recent calculations 
for HDO bv lWiesenfeld et al.l (|2011b ?) and for D2O bv lScribano et al.l (|2010[) . The HDO 
spectrum i s of particular inte rest as a probe of the depletion of molecules onto grain sur- 



faces (e.g., Parise et al.l 20051) . whereas D2O is only detectable under extreme conditions 



(JButner et al.lbOOTtlVastel et al.l |2010) 



Other notewort hy collision calculat ions from r ecent years are those for NH-He (jTobola et al. 



( Wiesenfeld et al 



2011a 



20111). CH+-He (iTurpinet al.l l2010l). SO2-H2 (|Cernicharo et al.ll201ll ). and NHD2-H 



In regions with strong radiation fields such as PDRs, collisions with electrons may 
contribute to the excitation of molecules. The calculation of collisional excitation proba- 
bilities wi th electrons is a sp ecialized field with its own approximations a nd pitfalls, as re- 
viewe d by lFaure et al. ( 2009f) . Examp les of recen t calculations are H j"-e~ ( Kokoouline et al. 
20101) . CS-e" (|Varambhia et al.l2010l). H9 0-e- (IZhang et al.l2009t) . SiO-e" (|Varambhia et al" 



20091 ) and the calculations bv lFaure etal] (|2007b[ ) for the excitation of HCN, HNC, DON 
and DNC in collision with electrons. 



3.4. H2O as an example 

The excitation of interstellar H2O is a long-standing problem in astrophysics, which 
remains to be solved today. At this moment, several calculations exist for its collision 
rates, which apply to different situations. At the low temperatures ( < 20 K) as occur 
for instance in dense pre-stellar cores and t he outer zones of proto planetary disks, the 
detailed quantum- mechanical calculations bv lDubernet et al.l ( 2006a ) are the best rates to 
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use. For higher temperatures (up to ^300 K), such as in mo deling protos t ellar en velopes 
and photon-dominated regions, the quasi-classical rates by iFaure et al.l ( 2007a ) are to 
be preferred, which at low temperatures join smoothly with the full quantum results. 
Finally at hi gh temperatures ( ^ 300 K) where vibrational excitation becomes important, 
the rates bv lFaure fc JosselinI (|2008[ ) should be used. 

The existence of one set of collision rates for H2O that appli es to all temperatures 



is obviously desir able, and work is ongoing towards this goal ((Dubernet et al.l bOOil 



Daniel et al.ll2010l ). Currently the quantum calculations for H2O at high temperature 



have only been published for collisions of ortho-H20 with para-H2 and para-H20 with 
ortho-H2, so that 50% of the rates is unavailable. Publication of the remaining rates 
woul d be very useful for the in odeling of H2O observations, in pa rticular with Hersc hel 
(e.g.. Ivan Dishoeck et al.ll2011 ). Use of the older H2 0-He rates bv iGreen et al.l ( 19931) is 
not recommended, as these numbers lie factors 3-4 below the corresponding rates with 
H2, even after scaling for the different masses. 

The H2O molecule is one of the rare cases where collision calculations have been verified 
in a laboratory experiment. While direct verification of state-to-state collision rates is 
limited to a few special circumstances (see § 13. 2p , pressure broadening experiments offer 
a way to measure the total inelastic collision rates to and from a given molecular energy 
level. Thus, measuring the derivative of the line width as a function of gas pressure is 
a way to test the accuracy of the potential energy surface that is also used to calculate 
c ollision rate s. 

Dick et al.l ( 20091 l2010r ) have measured the pressure broadening of the ground-state 



lines of ortho- and para-H20 in collisions with H2 gas as a function of temperature. They 
find that at temperatures above w80K, the slope of the pressure broadening curve as a 
function of tempe rature agrees well w ith theoretical predictions using the 5D potential 
energy surface bv lValiron et al.l ( 20081 ). and that the vertical level of the curve indicates 
an ortho-para ratio of w3 for the H2 gas. However, at temperatures below «80K, the 
experiments indicate a sharp drop of the pressure broadening coefficient, which is in 
appa rent contradiction with the theoretical predictions. 

As I Wiesenfeld fc Faurel ( 2010l ) point out, the disagreement between theory and experi- 
ment may only be apparent, because the interpretation of the experimental data depends 
critically on the approximatio n that the molecule s do not rotate during their collision. 
However, new experiments by iDrouin et al.l (|2011l ) show that this so-called 'impact ap- 
proximation' holds even in the case of collisions with H2 at low temperatures, and that 
instead, conversion of ortho- into para-H2 in the gas cell has caused the apparent decrease 
in the measured pressure broadening coefficients at T < 80 K. 

State to state experiments provide an even stro nger test o f a co llisional potential 



energy surface. In the case of the H2O-H2 system, I Yang et al.l (|2011l ) have achieved a 



confirmation of the quantum mechanical calculations using a crossed molecular beam 
experiment where the angular velocity distributions of fully specified rotationally excited 
final states are obtained using a technique called velocity map imaging. The experimental 
data agree with the theoretical predictions in impressive detail. 

3.5. What if collision data do not exist? 

Astronomers are often in the situation that collision data are unavailable for the molecular 
line that they have observed, so that they need workarounds. One option is to assume 
LTE and accept the limited accuracy of the model prediction, but often it is possible to 
use existing collision data with minor modifications. 

For isotopologues, adopting the collision data for the main species is a common strat- 



Radiative transfer and molecular data 9 

egy. This method works weU for species hke C'^^S, which have a similar structure as the 
main species, so that the potential energy surface of the collision is the same. The mass 
difference with the main species C'^^S is also small, so that the collision dynamics are 
similar. On the other hand, neglecting isotope effects is not an option for species like 
HDO, which lacks the ortho/para moieties of its main species H2O, and where the mass 
effect is also larger. 

Collision data for vibrationally excited states exist for very few molecules only. The 
coUisional excitation of vibrational modes is therefore often neglected, which is a good 
approximation if the radiative decay from these levels is fast, as is often the case. For 
low-lying vibrational states with small Einstein A coefficients, LTE models are the only 
option if collision data do not exist. 

For certain molecules, it is possible to borrow the collision data from other species 
with a similar structure. Such 'substituted' rates are generally good enough for non-LTE 
estimates of molecular column densities, but not good enough to constrain the kinetic 
temperature or the volume density of the parent gas. Cases where substitution works well 
include H2S and H2CS, where H2O and H2C O data are used resp ectively, after scaling the 
collision rates for the difference in mass (see Schoier et al.ll2005l for details) . Similarly, it 
should be possible to use scaled HE rates (JReese et al.ll2005l ) for HCl. The same approach 
is sometimes applied to cases like HOC+, where HCO+ rates are used, but the result is 
less accurate because the dipole moments of the two species differ significantly. 

Frequently, collision data only exist for low-excitation lines of a molecule, for a limited 
range of temperatures. To model observations of high-excitation lines of this molecule, or 
to model conditions outside the range for w hich collision data e xist, extrapolation of the 
existing data is the most common strategy. ISchoier et al.l (|2005l ) r uake recommendat ions 
on how to perform such extrapolations for small molecules, and iFaure et al.l (|2008l ) do 
so for larger molecules. 

Finally, cases exist where only collision data with He as partner are available. The 
common approach is to scale these rates to the different mass of the H2 molecule, which 
only takes care of the dynamics part of the interaction. Dedicated studies show that 
the potential energy surfaces of X-He and X-H2 sometimes differ substantially, espe- 
cially if collisions with ortho-H2 with its non-zero quadrupole moment are included (e.g., 
Cernicharo et al.ll201ir) . CoUisional excitation rates with H2 appear to be larger by fac- 
tors of 3-10 than expec ted by scaling He rates, a s indicated by several calcuations and 
also some experiments ( Mengel fc De Lucial 120001 ). 



4. Future developments 

In the coming years, the field of radiative transfer and molecular data will undergo 
rapid changes to keep up with the increasing data rates and higher precision of new 
generations of telescopes and instruments. Comparison of observations and models by 
hand will become impossible, and automated routines will appear which will search large 
pieces of parameter space for the best match to a set of observations. While such tools 
exist already at the individual level, the development of a parameter search engine for 
general use is actually a major undertaking. One such attempt is the CATS project (see 
the contribution by Schilke to these proceedings) 



A second, more ambitious project is ARTIST (jPadovani et al.ll2011al ). which besides 
a parameter search engine aims to develop a coupled radiative transfer program for 
dust and gas, where the thermal balance calculation is integrated with the radiative 
transfer. Like CATS, a link to solving chemical networks is foreseen, which is important 
for highly reactive species such as 0H+ and H2O+. Third, ARTIST aims to develop 
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continuum radiative transfer tools with polarization capability, which is of interest for 
ALMA observations which will be sensitive to polarization. 

Another future development is the merging and streamlining of the various databases. 
A few years ago, the Splataloguqj has appeared as a useful front-end to the JPL and 
CDMS spectroscopic databases, allowing for easy searching and switching between vari- 
ous units. One useful addition would be the option to create a state list and to calculate 
partition functions at any temperature. At some point, a compreh ensive database of spec - 
troscopic and collision information may appear as the VAMDC ( Dubernet et al.l 120101 ). 



While such meta-databases certainly make the access to molecular data easier for non- 
expert users, close ties between database developers and data sources are essential to 
ensure that the information in the database is up-to-date, and to ensure that the spec- 
troscopists, physicists, and chemists who produce the data get proper credit for their 
important work. 

Last but not least, the future will see new calculations of molecular collision rates. 
Commonly observed molecules for which no collision data currently exist include CH 
and C2H; for the latter, calculations are ongoing in Meudon and Le Havre. The col- 
lisional excitation of the hyperfine transitions of NH and ND is also being studied in 
Le Havre, while the Grenoble group is working on the H2O isotopologs. Other species 
for which collision calculations would be useful are 0H+ and H2 0"^, especially to model 



the o bservation of their emission lines in some extragalactic systems (jVan der Werf et al 



2010r) . For 'large' molecules (with more than Ri4) atoms, the quantum-mechanical calcu 



lation of collision potentials will likely remain impractical for the foreseeable future, and 
developing recipes or approximations for their coUisional excitation probabilities would 
be very useful for the interpretation of s pectral line surveys with Herschel and ground- 



based telescopes. See iFaure et al.l (|201lf ) for the example of HCOOCH3. Calculations 
of the excitation of atomic fine structure lines in collisions with H2 would be helpful 
to understand the p hysical structure of photon-dominated regions and 'warm dark gas' 



m 



(JLangeret al.l 120101) : for O, C, N+ and C+, only fairly old («1990) data exist, while 



for S, no calculations exist at all. Finally, inclusion of electronic and vibrational energy 
levels and radiative transitions into molecular databases will be crucial to quantify the 
radiative excitation rotational lines through infrared, optical and ultraviolet pumping. 
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Discussion 

HORACEK: My group has measured the rovibrational and electronic spectra of HF and 
HCl, which may be useful in include in the astrophysical data bases. 

Van der Tak: Yes, those data would be very useful. 



